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ABSTRACT 


The nearby X-ray binary X Per (HD 24534) provides a useful beacon with 
which to measure elemental abundances in the local ISM. We examine absorption 
features of 0, Mg, and Si along this line of sight using spectra from the Chandra 
Observatory’s LETG/ACIS-S and XMM-Newton’s RGS instruments. In general, 
we find that the abundances and their ratios are similar to those of young F and G 
stars and the most recent solar values. We compare our results with abundances 
required by dust grain models. 


Subject headings: ISM: abundances 



1. Introduction 


An integral part of understanding the interstellar medium (ISM) and building realistic 
du3t grain models is determining the abundances of the elements in the dust phase. The 
most common method is to measure this indirectly, that is, to measure the gas-phase 
abundance and subtract this from the total (gas+dust) abundance. The gas-phase 
component of the local ISM is fairly well-known, having been determined by UV /optical 
absorption studies by many workers (Cardelli et al. 1996; Sofia et al. 1997; Cardelli & 
Meyer 1997; Sofia et al. 1998; Jensen et al. 2005; Gnacinski & Krogulec 2006) and appears 
to be relatively independent of the line of sight used. However, the total local ISM elemental 
abundances are not well known, and proxies must be used instead. 

Historically, the Sun has been used as the standard for ISM abundances. However, 
while most solar abundances have remained more or less constant over the years, those 
of the main dust-forming elements have been revised downward considerably (Anders & 
Grevesse 1989; Grevesse & Sauval 1998; Asplund et al. 2005; Asplund et al. 2009). The 
recent revisions have brought solar O into better agreement with inferred ISM abundances 
of Wilms et al. (2000) and Sofia & Meyer’s (2001) young F and G stars. We note, though, 
that large scatter remains and there is no particularly persuasive reason why the 4.5 Gyr 
Sun should be an accurate representation of the ISM today. 

B stars and metal-rich F and G stars have also been suggested as proxies for the ISM 
■but these have problems of their own. Hempel k Holweger (2003) showed that B stars 
cannot be considered an accurate indicator of ISM abundances, as element stratification 
due to diffusion appears common in these stars. Thus the abundances measured at their 
surfaces do not necessarily reflect the abundances of the clouds from which they formed, 
though the latest revision of B star abundances seem to ease differences between them and 
solar values somewhat (Nieva k. Simon-Diaz 2011). Regarding F and G stars, Edvardsson 
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et al.’s (1993) study of stars in the solar neighborhood showed that there is large scatter in 
the their age-metallicity relation; thus,- metal-rich F and G stars are not necessarily young, 
and cannot be used as a gauge of today’s ISM abundances. Sofia & Meyer (2001), using 
Edvardsson et al.’s (1993) sample, calculated the average abundances for F and G stars less 
than 2 Gyr old. They found values for Mg and Si that are similar to the most recent solar 
values, but less 0 than solar (Asplund et al. 2009) Some have suggested using abundances 
measured in H II regions, but this has been rejected due to large systematic errors from 
temperature fluctuations in the nebula (Mathis 1995; Kingdon & Ferland 1995). Further 
complicating the situation, and highlighting the need for spectroscopy on nearby objects, 
is the correlation between Galactocentric distance and abundances (Shaver et al. 1983; 
Gummer sbach et al. 1998; Chiappini et al. 2001), as regions closer to the Galactic center 
are more metal-rich than outer regions. 

A direct measurement of the local ISM is the best way to determine exactly how much 
material remains for local dust grains. X-ray spectroscopy is particularly useful for this, as 
the elements which make up the grains have their K- and L- shell absorption edges in this 
band. The edge depths are dependent on each element’s total column density (gas + dust 
phases), while fine structure near the edges can determine the chemical state. Schattenburg 
&; Canizares (1986) were the first to measure X-ray absorption edges; they used the Einstein 
Observatory to study the K edges of O and Ne toward the Crab Nebula. Since then, much 
X-ray spectroscopy has been done, usually toward very bright (and distant) low mass X-ray 
binaries. Paerels et al. (2001) used the Chandra X-ray Observatory’s LETG to measure 
the column densities of O, Ne, and Fe toward 4U0614+091, at a distance of about 1.5 
kpc (Machin et al. 1990). Takei et al (2002) examined the O edge and its fine structure 
along the line of sight toward Cyg X-2, and Ueda et al. (2005) determined the abundances 
of O, Mg, Si, S, and Fe toward GX 13+1, GX 5-1, and GX 340+0. Juett et al. (2004, 
2006) examined the K edge and fine structure of O, Ne, and Fe for several x-ray binaries. 
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Constantini et al. (2010) used XMM-Newton, INTEGRAL, and Chandra to study the ISM 
toward 4U1820-30, focusing on fitting fine structure near the Fe L and O K edges with dust 
constituent spectra derived from experiment. However, these studies focus on long lines of 
sight. For instance, 4U1820-30 is located in the globular cluster NGC 6624, between 6.8-9. 6 
kpc away (Guver et al. 2010). Cyg X-2 is 1.6 kpc away (Goranskij &, Lyutyj 1988), and 
GX 13+1, GX 5-1, and GX 340+0 are near the Galactic center, 7, 9.2, and 11.0 kpc away, 
respectively (Jonker et al. 2000). 4U0614+091 is about 3.2 kpc away (Kuulkers et al. 2010). 

In contrast, Cunningham et al. (2004) used the Chandra LETG to examine the line of 
sight toward X Per, only ~ 400 pc away. While they were able to measure the O column 
density, the observation unfortunately had low signal-to-noise, leading to an uncertainty 
of about 50% in their measurement and an inability to impose meaningful constraints on 
grain models. Thus, to better assess the elelmental budget in the nearby ISM, in particular 
that of silicate-forming dust, we have obtained and analyzed high-resolution X-ray spectra 
of X Per from both Chandra and XMM-Newton. 


2. The X Per System 

X Persei (HD 24534, BD+30 591) is a Be/X-ray binary is composed of a neutron star 
and OB companion of likely spectral type 09.5III to BOV (Liu et al. 2006; Raguzova 2007). 
It is variable, with my = 6.1 - 6.8. The optical variations are believed to be due to the 
formation and dissipation of an equatorial disk (Telting et al. 1998), similar to other such 
Be binary systems. The neutron star has a rotation period of about 835 s, one of the longest 
observed for a neutron star (Haberl et al. 1998). Distance estimates from spectroscopic 
parallax have placed it as far away as 1.3 kpc (e.g. Fabregat et al. 1992) but more recent 
work indicates it is closer, at ~ 0.7 - 1 kpc (Lyubimkov et al. 1997; Roche et al. 1997; 
Telting et al. 1998). Hipparcos parallax places it at 427lgg 2 3 pc, in agreement with Blaauw’s 
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(1952) assessment that it is in the £ Per OB association (d = 350 pc; Johnson 1959). 

The ISM toward X Per has also been studied extensively, and it happens that this sight 
line is an excellent representative for the diffuse ISM at large, as seen through the ratio of 
total- to-selective extinction Ry (=Ay/(E(B-V)). Ry is a general indicator of grain size and 
local ambient radiation field; it is the sole variable upon which differences in UV extinction 
curves depend (Cardelli et al. 1989). It ranges from ~2.5 (diffuse sight line) to ~5.5 (dense 
sight line), with the “standard diffuse ISM” having Ry= 3.05±0.15 (Whittet 2003); most 
grain models assume this standard value. It is identical to Ry for the X Per line of sight, 
3.06 (Valencic & Smith 2008), making X Per ideal for studying the local ISM. 


3. Observation and Data Reduction 

We used two observations, one from Chandra’s LETG/ACIS-S and the other from 
XMM-Newton’s RGS. The Chandra observation (Obs ID 12447) had an exposure time of 131 
ks. We reprocessed the data with CIAO 4.3 (CALDB 4.4.1) according to standard methods 
for grating spectra as described in the Science Threads 1 using the task “chandra_repro” to 
generate a new bad pixel file, event file, and PHA files. Before continuing, we examined 
the light curve and found that it did not show any flaring during this time, and only 
slight variability, with an average count rate of 0.96 ± 0.10 count/sec. We then created 
and applied the response and ancillary files for the +1 and -1 orders, which we coadded 
and grouped to ensure a minimum of 15 counts per bin (Cunningham et al. 2005), which 
led to a binsize of 0.007 keV over the spectrum. The XMM-Newton observation (Obs ID 
0600980101) had an exposure time of 126 ks. We reprocessed the RGS data with SAS 11, 


1 http: / / cxc.cfa.harvard.edu/ciao / 
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according to standard procedures described in the “ABC Guide” 2 . After reprocessing, we 
examined the light curve, which showed flaring at the start and end of the observation; once 
we excised those events, the effective exposure time was '-*'100 ks. We made and applied 
response files for the first order spectra and rebinned to ensure a minimum of 15 counts per 
bin, which led to a binsize of 0.0013 keV over the spectrum. We did not make ancillary 
files, as these are not needed for RGS data. We also did not combine the spectra from the 
two RGS instruments. In the analysis, we examined the LETG data first, as it typically has 
a higher effective area than RGS. 

Unfortunately, the signal-to-noise was too low to reliably detect fine structure. 

The edges were measured as described below. We used the cross sections of Verner & 
Yakovlev (1995) and the total H column density from Cunningham et al. (2004), Nr = 
2.19xl0 21 cm~ 2 , to calculate the number of atoms per million H nuclei. The uncertainties 
given with our results are the 1 a statistical errors. Table 1 lists our results and the amount 
of material in the gas and dust phases, where the gas phase abundances are from the 
literature. Fig. 4 compares our results to those of other workers. 


3.1. Oxygen 

We began by examining the LETG data following the method described by Juett et al. 
(2004) , fitting the continuum around the edge with an absorbed power law, and freezing the 
continuum while fitting the edge explicitly. This was done over a range of energies (Takei et 
al. 2002). We used the tbvarabs model (Wilms et al. 2000), holding Nh = 2.19xl0 21 cm -2 
and O abundance at 0, and measuring the actual O absorption with an additional edge 
model. While the spectrum had suggestions of fine structure, in particular for the atomic O 


2 htt p : / /heasarc . nasa. gov / do cs /xmm/ ab c/ 
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absorption line at 0.527 keV, the signal-to-noise was too low to fit the feature. Therefore, 
we restricted ourselves to a simple edge fit and did not include the 0.525-0.548 keV region 
in the continuum fit (Takei et al. 2002; Cunningham et al. 2005). This also minimized 
confusion with instrumental effects. We tried to avoid absorption features from other 
elements as much as possible, specifically from the Fe L II and III edges at 0.71 and 0.72 
keV, respectively, and Ne at 0.87 keV. While no absorption from Fe L edges was visible, 
we nonetheless tested the dependence of the O optical depth on the Fe abundance by 
assuming Fe abundances of 0.5, 1.0, and 2.0 x the values of Wilms et al. (2000) (Takei et 
al. 2002) and following the same fitting procedure described above. This led to a maximum 
difference in optical depth of about 4%, or 30% of the statistical error. We therefore kept 
the Fe abundance fixed at the value from Wilms et al. (2000). We found that fits with 
a minimum energy of 0.44-0.46 keV and a maximum energy of 0.69-0.73 keV consistently 
produced edge energies and optical depths that were within 1 a of each other, with edge 
energies of 0.537±0.010 keV to 0.541±0.010 keV, and optical depths between 0.54±0.03 
and 0.57±0.03. Values of the reduced xl ranged from 0.68 to 0.72, with the best values 
corresponding to optical depths of 0.54-0.56 and an edge energy of 0.538±0.010 keV. 

We also examined the edge in the RGS spectrum, following the same procedure. We 
used only RGSl, as CCD4 in RGS2, which covers 0.514-0.620 keV, was lost early in the 
mission. Over the same energy ranges, we found the edge energies of 0.538±0.010 keV to 
0.540 ±0.010 keV. The optical depths arid xl values were in the same ranges as for LETG, 
and the best fits corresponded to optical depths of 0.54-0.55±0.03. Using an optical depth 
of 0.54±0.02, we find N 0 = (1.09 ± 0.04) x 10 18 cm- 2 , and N 0 /N H = (4.95 ± 0.18) x 10~ 4 . 
Our value for the edge energy is lower than expected for O 2 and atomic oxygen (0.543 and 
0.544 keV, respectively; Sevier 1979; Stolte et al. 1997) and in agreement with that of metal 
oxides (0.532~0.540 keV; Sevier 1979; Nakai et al. 1987) and molecules containing O-C 
single bonds (0.538-0.542 keV; Yang et al. 1987). Our value for No is slightly higher than 
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those found by Cunningham et al. (2004), though we still agree within 3 a; they found ~ 
7-8^4 ;q x 10 17 cm -2 . The fits from 0.45-0.70 keV are shown in Fig. 1. 


3.2. Magnesium 

We began by examining the LETG data, which does not have internal Mg features, 
unlike the RGS. Laboratory spectra of MgO show structure between 1.31-1.33 keV (i.e. 
Luches et al. 2004, Limpijumnong et al. 2011), with the K edge near 1.31 keV. No such 
struct me was visible, but we nonetheless excluded this region from our fits to minimize 
any effect it might have. We followed the same method as for O: we fitted the continuum 
around the edge, exluding the region from 1.31-1.33 keV, over a range of minimum and 
maximum energies. We again let the continuum be an absorbed power law, with Nh — 
2.19xl0 21 cm -2 in tbvarabs and the Mg abundance set to 0. For each energy range, we froze 
the continuum and tried to fit the edge, but the fits consistently could not find it. We then 
included the 1.31-1.33 keV range and allowed the continuum parameters to float; again, 
the fits could not consistently find the edge. We then froze the edge energy at 1.307 keV 
(Ueda et al. 2005) and re-fitted the spectrum. We found that for minimum energies from 
1.25-1.27 keV and maximum energies from 1.36-1.38 keV, the optical depth ranged from 
0.016-0.024±0.008, with the best fits (x^=0.94) consistently producing an optical depth of 
0.020±0.007. 

We then analyzed the RGS data. The RGS2 has a chip edge at 1.24 keV, near the Mg 
K edge, so we restricted ourselves to the RGSl to minimize any effect this might have on 
the fit. We again let the continuum be an absorbed power law, with Nh = 2.19x 10 21 cm -2 
in tbvarabs and set the Mg abundance to 0. We again began by fitting the continuum 
around the edge, exluding the region from 1.31-1.33 keV, over a range of minimum and 
maximum energies. For each energy range, we then froze the continuum and tried to fit the 
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edge, but this consistently yielded optical depths smaller than the uncertainies. Despite the 
possibility of introducing instrumental effects, we then included the 1.31-1.33 keV region in 
the fits and allowed the continuum parameters to float. For the same range of minimum 
and maximum energies, we consistently found an edge energy that was within 1 a of the 
LETG result, 1.309±0.007, and the best fits (x 2 =1.19) corresponded to an optical depth of 
0.G19±0.013. For an optical depth of 0.020±0.007, we find NM g = (9.86 ± 3.45) x 10 16 cm -2 
and NMg/Nn = (4.50 ± 1.58) x 10 -5 . The fits from 1.25-1.38 keV are shown in Fig. 2. 


3.3. Silicon 

We again fitted the LETG spectrum by treating the continuum as an absorbed power 
lav with Nh at 2.19xl0 21 cm~ 2 in tbvarabs and setting the Si abundance to 0 in tbvarabs. We 
excluded the region 1.84-1.87 keV in an effort to minimize any confusion with instrumental 
effects (Ueda et al. 2005), and allowed all parameters to float. We found that fits with a 
minimum energy 1.75-1.78 keV and maximum energy 1.93-1.95 keV produced the best fits, 
with xl between 0.97-1.03 and values for the optical depth of 0.007±0.006 to 0.010±0.006 
and edge energy of 1.846±0.015 to 1.851±0.015. The best values of xl corresponded to an 
optical depth of 0.010±0.006. The fit from 1.75-1.95 keV is shown in Fig. 3. 

We again tried to verify this measurement with the RGS data using the same energy 
ranges and method. We fitted the R1 and R2 spectra simultaneously, and found no excess 
absorption at the Si edge. However, the R1 and R2 have bad pixels near the Si edge, which 
could lead to underestimates of the continuum. To see if these were having an effect on the 
measurement, we then fitted the R1 and R2 spectra seperately. Again, no excess absorption 
was detected. 


We then attempted to verify the origin of the additional absorption (ISM or instrument) 
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in the Chandra spectrum by examining an unabsorbed line of sight. We examined a 
Chandra LETG/ACIS-S spectrum of Mkn 421 (Obs ID 5171), a line of sight known to 
have very little interstellar material. We used the HEASARC’s Galactic HI column density 
calculator to estimate Nh within a 1° radius of Mkn 421. The Galactic HI survey of Dickey 
& Lockman (1990) indicates Nh ~ 1.4xl0 20 cm -2 , and the Leiden/Argentine/Bonn Survey 
(Kalberla et al. 2005) shows Nh ~ 1.9xl0 20 cm" 2 , both of which are lower than that of X 
Per’s line of sight by an order of magnitude. We reprocessed the Mkn 421 data, extracted 
the spectrum, and examined it about the Si edge using the same methods and over the 
same energy ranges as the X Per spectrum. No excess absorption was found, so despite 
the smallness of the measurement for the X Per line of sight, we believe it represents 
the contribution from the ISM. Thus, for an optical depth of 0.010±0.006, we find N Si = 
(7.37±4.42)xl0 16 cm" 2 and N Si /N H = (3.37 ±2.02) x 10" 4 5 . 


4. Discussion 

In general, our results tend to be similar to those found by other studies of the local ISM 
and the most recent proto-solar nebula estimates, as determined by the diffusion-corrected 
photospheric values of Asplund et al. (2009), which has brought solar 0 values closer to 
those of young F and G stars (Sofia &: Meyer 2001). Fig. 4 shows our results and those of 
other workers for comparison. Abundances of O show a fair amount of scatter, from 251 

atoms/10 6 H in meteorites (Lodders et al. 2009) to 851 atoms/10 6 H in the solar value 
from Anders & Grevesse (1989). However, recent studies show 0 abundances between these 
values. A re-analysis of 13 B stars in Ori OBI by Nieva & Simon- Dfaz (2011) indicated 
there are 589 0 atoms/10 6 H, though another study by Simon-Dfaz (2010) places the O 
abundance for the same 13 stars at 537 atoms/ 10 6 H. Our value of 495 atoms/ 10 6 H is 
nearly identical to that of Wilms et al. (2000) and agrees within 3 a of Sofia & Meyer’s F 
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and G stars, Simon-Dfaz’s (2010) B stars, and Asplund et al.’s (2009) solar nebula estimate. 

There is not nearly so much scatter in Mg estimates, which tend to range from about 
25 to 40 atoms/10 6 H. While our result for Mg is high (45.0±15.8), it is nonetheless 
similar to that of Asplund et al. (2009) and Sofia & Meyer (2001) (43.7l3;g and 42.7+17.2, 
respectively). Estimates for Si also tend to cover a narrow range, about 20 to 40 atoms/ 10 6 
H. As with Mg, our result is similar to that found by Asplund et al. (2009) and Sofia & 
Meyer (2001), as well as Simon- Diaz (2010). 

Seeing the agreement between our abundance determinations and those of other 
workers, in particular that of Asplund et al. (2009), Sofia & Meyer (2001), and Simon-Diaz 
(2010), it is not surprising that our abundance ratios are similar to theirs, too. Fig. 5 
compares our ratios to those of the same workers from Fig. 4. 

We also compared our results to those made toward other X-ray binaries. Juett et 
al. (2004), measured the O column density and examined structure near the O edge in 
seven lines of sight toward bright X-ray binaries. In order to estimate the total Nh, we 
made use of Takei et al.’s (2002) finding that for the diffuse ISM (Nhi = 2.17 x 10 21 cm -2 ) 
toward Cyg X-2, Nh 2 ~ 5% Nhi and Nhii ~ 32% Nhi. We again used the Galactic HI 
column density calculator and the Leiden/Argentine/Bonn Survey (Kalberla et al. 2005) 
to find N H i within 1° for each of the sources, then found the total H column density. This 
resulted in O abundances from 378-916 atoms/10 6 H. Unsurprisingly, our abundances are 
notably lower than those along the sight lines toward the Galactic Center (Ueda et al. 
2005). Using Ueda et al.’s (2005) total Nh values for GX 13+1, GX 5-1, and GX 340+0, 
(2.1Z.I I x 10 22 , 2.81^8 x 10 22 , an d 4.8±2:j) x 10 22 cm -2 , respectively), we found that 0 ranged 
from 833-857 atoms/10 6 H, Mg ranged from 75-86 atoms/10 6 H, and Si ranged from 60-75 
atoms/10 6 H. In the case of GX 5-1, Smith et al. (2006) found that Nh could be as high 
as 3.3 xlO 22 cm -2 , depending on the brightness of a CO cloud along the sight line. This 
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higher Nh, while lowering the abundances commensurately, did not lower them, sufficiently 
to bring them into agreement with values from the local ISM: 727, 72, and 55 atoms/10 6 
H for O, Mg, and Si, respectively. Our 0 abundance is also substantially lower than that 
which Costantini et al. (2012) found for the line of sight toward 4U1830-20, 755 atoms/10 6 
H. 

Using our total abundances, and those from the gas phase from the literature, we can 
find how much is depleted into dust. These are shown in Table 1. Only about 38% of the 
total O for this line of sight is in the dust phase, while Mg and Si are almost completely 
depleted (97%). This is in accordance with what has been found for other diffuse lines 
of sight (Savage & Sembach 1996; Whittet 2003). Snow & Witt (1996) examined the 
abundance requirements of nine grain models; these are fisted in their Table 3. In a similar 
manner, we examine more recent models in our Table 2, and include their findings. We 
also include the dust phase abundances found in this work for ease of comparison. For 
models where only Si is fisted, we assumed that the O requirements were four times that 
of Si (SiO,i) (Snow & Witt 1996; Cunningham et al. 2004) and for models where Mg was 
not listed and the type of silicate was not specified, we assumed that the model needed an 
amount comparable to Si (Snow & Witt 1996). Abundance requirements found in this way 
are in parentheses. 

For this fine of sight, there appears to be sufficient O for most models to within 1 
a. Unfortunately, the large uncertainties in our Mg and Si measurements do not allow us 
to meaningfully constrain these quantities. We can only note that there appears to be 
sufficient Mg for most models, and less Si appears available than needed by the models of 
Weingartner & Draine (2001) and Iati et al. (2008). 
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5. Conclusions 

We have obtained and analyzed spectra of the short line of sight toward X Per using the 
Chandra and XMM-Newton X-ray observatories. X-ray spectroscopy shows the abundances 
of silicate dust-forming elements are largely similar to the latest estimates from classical 
proxies. The 0 abundance is substantially higher than that determined from meteorites. 
While the uncertainties in the Mg and Si measurements are large, our findings suggest 
that the abundances of these elements are more similar to solar, young F and G star, and 
the latest B star estimates than previously inferred ISM values. Higher signal-to-noise 
observations are needed to reduce the large uncertainties in the Mg and Si measurements 
and place meaningful constraints on dust grain models. 

Facilities: CXO (ACIS), XMM-Newton (RGS) 
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Fig. 2. — Best fits for Mg for LETG (left) and RGS (right). The gap in the RGS spectrum 
at 1.36 keV is from a bad pixel. The location of the edge is indicated by the dotted line. 




Normalized 


20 



Fig. 3. — Best fit for Si for LETG. The location of the edge is indicated by the dotted line. 
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Fig. 4. — Number of 0, Mg, and Si atoms per 10 6 H nuclei. Results from other studies 
are shown for comparison. AG89=Anders & Grevesse 1989; A09= Asplund et al. 2009; 
SM01=Sofia & Meyer 2001; SD10=Simon-Diaz 2010; W00=Wilms et al. 2000; L09=Lodders 


et al. 2009. 
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Fig. 5. — Ratio of 0, Mg, and Si abundances. Results from other studies are shown for 
comparison. AG89=Anders & Grevesse 1989; A09= Asplund et al. 2009; SM01=Sofia & 
Meyer 2001; SD10=Simon-Diaz 2010; W00=Wilms et al. 2000; L09=Lodders et al. 2009. 
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Table 1: Gas and Dust Phase Abundances of 0, Mg, and Si Relative to 10 6 H Nuclei 



Total 

Gas phase 0 

Dust phase 

0 

495.4±18.3 

308 . 7 ^ 9.8 

186 . 71^6 

Mg 

45.0±15.8 

1.24±0.09 

43.8±15.8 

Si 

33.65±20.19 

0.92±0.55 

32.73±20.19 


°0: Jensen et al. 2005; Mg and Si: Gnacinski & Krogulec 2006 
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Table 2: Comparing Grain Model Requirements per 10 6 H Nuclei and the Dust Phase Abun- 
dances 


Reference 
Snow & Witt 1996 
Weingartner & Draine 2001 
Zubko et al. 2004 
Vaidya et al. 2007 
Iati et al. 2008 
Cecchi-Pestellini et al. 2010 


O Mg Si 

124-144 20-36 27-40 

182-193 45-48 45-48 

114-176 28-34 28-34 

(96-112) (24-28) 24-28 

(180) (45) 45 

(136) (34) 34 

187 44 


This work 
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